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Recent measurements of temperature fluctuations in the cosmic microwave background 
(CMB) indicate that the Universe is flat and that large-scale structure grew via gravi- 
tational infall from primordial adiabatic perturbations. Both of these observations seem 
to indicate that we are on the right track with inflation. But what is the new physics 
responsible for inflation? This question can be answered with observations of the polar- 
ization of the CMB. Inflation predicts robustly the existence of a stochastic background 
of cosmological gravitational waves with an amplitude proportional to the square of the 
energy scale of inflation. This gravitational-wave background induces a unique signature 
in the polarization of the CMB. If inflation took place at an energy scale much smaller 
than that of grand unification, then the signal will be too small to be detectable. How- 
ever, if inflation had something to do with grand unification or Planck-scale physics, 
then the signal is conceivably detectable in the optimistic case by the Planck satellite, or 
if not, then by a dedicated post-Planck CMB polarization experiment. Realistic devel- 
opments in detector technology as well as a proper scan strategy could produce such a 
post-Planck experiment that would improve on Planck's sensitivity to the gravitational- 
wave background by several orders of magnitude in a decade timescale. 

1. What Have We Learned from the Cosmic Microwave Background? 

The past year has seen spectacular advances in measurements of temperature fluc- 
tuations in the cosmic microwave background that have led to major 
advances in our ability to characterize the largest-scale structure of the Universe, 
the origin of density perturbations, and the early Universe. The primary aim of 
these experiments has been to determine the power spectrum, Ci, of the CMB as a 
function of multipole moment I. Given a map of the temperature T(n) in each di- 
rection n on the sky, the power spectrum can be obtained by expanding in spherical 
harmonics, 

a tm = J dhY em (h)T(h), (1) 
and then squaring and summing the coefficients, 

^ = ^TlEl a ^l 2 - ( 2 ) 
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If the map covers a patch of the sky that is small enough to be approximated as a 
flat surface, the power spectrum can be written in terms of Fourier coefficients: 

T ( = J dhe- lle ~T(h), (3) 

and then 

C t-(\ T d%\=t> ( 4 ) 

where the average is taken over all Fourier coefficients I that have amplitude I. 
Thus, each multipole moment Ce measures, roughly speaking, the rms temperature 
fluctuation between two points separated by an angle 9 ~ (f/200) _1 degrees on the 
sky. 

Recent experiments have sought to determine the power spectrum in the range 
50 < £ < 1000, as structure-formation theories predict a series of bumps in this 
regime that arise as consequences of oscillations in the baryon-photon fluid in the 
era before recombination (as indicated by the curves in Fig. |l|). The rich structure 
in these peaks, which can be characterized, e.g., by the precise heights and widths 
of the peaks, their locations in £, and the heights of the troughs between the peaks, 
depends in detail on the values of several classical cosmological parameters, such 
as the baryon density f2f, (in units of the critical density), Hubble constant h (in 
units of 100 km/sec/Mpc), matter density f2 m , and cosmological constant on 
structure-formation parameters such as the amplitude and spectral index of pri- 
mordial perturbationsB'El; and on the character of primordial perturbations (e.g., 
adiabatic, isocurvature, or topological-defects products). In particular, the location 
of the first peak depends primarily on the geometry of the Universe (parameter- 
ized by the total density O to t)) and only secondarily on the other cosmological 
parameters!! If the Universe is flat, the first peak is expected to occur at I ~ 200, 
while if the Universe has a matter density fl rn ~ 0.3 (as dynamical measurements 
indicate) but is open (no cosmological constant), then the first peak should be at 
I ~ 500. 

Experiments that measure the power spectrum in the regime 50 < I < 1000 
require high sensitivity to detect the CMB temperature variations of roughly one 
part in 100,000, and they require subdegree angular resolution. Within the past 
year, the first high-signal-to-noise high- angular- resolution maps of the CMB have 
been published bythe BOOMERanGfl and MAXIMA0 collaborations. The results 
of a joint analysis^ of the data from both BOOMERANG and MAXIMA are shown 
in Fig. [l]. The data show a peak at i ~ 200 which provides very strong evidence 
that the Universe is flat (earlier measurements by the TOCO collaboration!!! also 
indicated a first peak at I ~ 200, but with lower signal-to-noise). The peak structure 
is also very consistent with growth of large-scale structure from a nearly scale- 
invariant spectrum of primordial adiabatic perturbations and very inconsistent with 
isocurvature or topological-defect alternatives. The peak structure indicated in 
Fig. |l| is also beginning to provide valuable information about the values of other 
cosmological parametersou. 
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Fig. 1. Data points on the CMB temperature power spectrum obtained individually by 
BOOMERanG and MAXIMA, as well as the points obtained from a joint analysis of the two 
data sets. The data indicate unequivocally a peak at £ ~ 200 and are beginning to show the 
outline of a second peak at £ ~ 500. From Jaffe et alo 
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2. Inflation and Gravitational Waves 

The flatness of the Univejseand adiabatic perturbations suggest that we are on the 
right track with infLatior£3lld'E3, a period of accelerated expansion in the very early 
Universe driven by the vacuum energy associated with some new ultra-high-energy 
physics. In order to solve the horizon problem for which it was initially proposed, 
inflation predicts that the Universe is flat. Moreover, shortly after inflation was 
proposed, it was realized that vacuum fluctuations in the inflaton (the scalar field 
responsible for inflation') would produce a nearly-scale-invariant spectrum of adia- 
batic perturbationsEi00'0O. With the advent of these CMB tests, inflation has 
now had several opportunities to fail empirically, but it has not. Conservatively, 
these successes are at least suggestive and warrant further tests of inflation. 

Perhaps the most promising avenue toward further tests of inflation is the 
gravitational-wave background. In addition to predicting a flat Universe with 
adiabatic perturbations, inflation also predicts that quantum fluctuations in the 
spacetime metric during inflation would give rise to a_atochastic gravitational- 
wave background with a nearly-scale-invariant spectrumlla. Quantum fluctuations 
in the spacetime metric can only be affected by gravitational effects which are 
quantified completely during inflation by the expansion rate H ln f[. This is related 
through the Friedmann equation to the vacuum-energy density V during inflation, 
H? nf[ = 87rVy(3TOp[), where mpi is the Planck mass (mpj 2 = G, Newton's constant, 
in particle-physics units h = c = 1). Thus, the amplitude of the gravitational- wave 
background is fixed entirely by the vacuum-energy density during inflation, which 
itself should be proportional to the fourth power of the energy scale E- ln s of the new 
physics responsible for inflation. The spectrum of gravitational waves depends on 
the particular inflationary model, but in most models (and certainly in the simplest 
inflationary models), it is likely to be very close to scale invariant. 

These gravitational waves will produce temperature fluctuations at large angles 
(I < 1 100) in the CMB (as shown in Fig. ^). The amplitude of their contribution 
to the CMB temperature quadrupole can be writtenH 

T = 6C 2 TT ' tcnsor = 9.2-^, (5) 
m pl 

(where "tensor" refers to gravitational waves, as they are tensor perturbations to 
the spacetime metric and "TT" refers to the temperature quadrupole). Since the 
quadrupole measured by COBE, Cj T — (1.0 ± 0.1) x 10~ 10 , is most generally due 
to some combination of density perturbations and gravitational waves, we already 
have animportant constraint to the energy scale of inflation: U 1 / 4 < 2 x 10 16 

GevSii 

3. Gravitational Waves and Polarization 

But how can we go further? One might think that improved temperature maps 
could be used to measure the power spectrum well enough to distinguish the rela- 
tive contributions of the gravitational- wave and density-perturbation power spectra 
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Fig. 2. Temperature and polarization power spectra from density perturbations (dotted curves) 
and gravitational waves (solid curves). The absence of a dotted curve for the CC (lower right- 
hand) panel is due to the fact that density perturbations do not produce a curl component. The 
solid curves show predictions for a model in which there is no reionization. More realistically 
some fraction r ~ 0.1 of the CMB photons will have re-scattered from reionized gas, and this will 
generate additional polarization power at large angles, as indicated by the dashed curve in the CC 
panel. 
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indicated in the upper left panel of Fig. ^. However, the precision with which the 
power spectrum can be measured is limited even in an ideal experiment by cosmic 
variance, the sample variance due to the fact that we have only 2i + 1 independent 
modes with which to measure each Ce- 

Instead, progress can be made with the polarization of the CMB. In addition to 
producing temperature fluctuations, both gravitational waves and density pertur- 
bations will produce linear polarization in the CMB, and the polarization patterns 
produced by each differ. This can be quantified with a harmonic decomposition of 
the polarization field. The linear-polarization state of the CMB in a direction n can 
be described by a symmetric trace-free 2x2 tensor, 

1 / Q(n) -U(n)sm6 \ 
/ afc (n) - - ^ _ u{h)sin9 _ Q{h)sin 2 e J , (bj 

where the subscripts ab are tensor indices, and Q(n) and U(h) are the Stokes 
parameters. Just as the temperature map can be ex panded . in terms of spherical 
harmonics, the polarization tensor can be expandedJ2jr 2 tr 3 H 2 'i 

^f—^ = X/ [ a (in)^>n)afc(n) + a (lL)^(?m)ab(n) > CO 
Im 

in terms of tensor spherical harmonics, Y{f m ^ ab and Y^ m \ ab ' It is well known that a 
vector field can be decomposed into a curl and a curl-free (gradient) part. Similarly, 
a 2 x 2 symmetric traceless tensor field can be decomposed into a tensor analogue 
of a curl and a gradient part; the Y^f m ^ ab and Yy m \ ab form a complete orthonormal 
basis for the "gradient" (i.e., curl- free) and "curl" components of the tensor field, 
respectively. * Lengthy but digestible expressions for the Y^f m ^ ab and Y9 m \ ab are 
given in terms of derivatives of spherical harmonics and also in terms of Legendre 
functions in Kamionkowski et aLc 2 ] The mode amplitudes in Eq. (^) are given by 

4m) = Y J dh V ab (n) Y^ b * (n) , 

«(L) = YjdhV ab (h)Y^*(h), (8) 

which can be derived from the orthonormality properties of these tensor harmonicsEl. 
Thus, given a polarization map 'P (,(n), the G and C components can be isolated by 
first carrying out the transformations in Eq. (^|) to the a^ m ) and o^lmV anc ^ ^ nen 
summing over the first term on the right-hand side of Eq. (|?]) to get the G compo- 
nent and over the second term to get the C component. (In practice, a full likelihood 
formalism would be used to determine the spectra in the presence of anisotropic, 
correlated noise, astrophysical foregrounds, and incomplete sky coverage.) 

The two-point statistics of the combined temperature/polarization (T/P) map 
are specified completely by the six power spectra C^ x for X, X' = {T, G, C}. Parity 

■fOur G and C are sometimes referred to as the "scalar" and "pseiidn-scalar" componentsEl , re- 
spectively, or with slightly different normalization as E and B modesEj (although these should not 
be confused with the radiation's electric- and magnetic- field vectors). 
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invariance demands that Cj c — CfSL = (unless the physics that gives rise to 
CMB fluctuations is parity breakingc3E3) . Therefore, the statistics of the CMB 
temperature-polarization map are completely specified by the four sets of moments: 

^TT WTG ^ GG and r cc 

Both density perturbations and gravitational waves will produce a gradient 
component in thp_p0larization. However, only gravitational waves will produce 
a curl componentO>E3 (but see below) . The curl component thus provides a model- 
independent probe of the gravitational-wave background, and it is thus the CMB 
polarization component that we focus on here. 

4. Detectability of the Curl Component 

If our goal is to detect the polarization signature of gravitational waves, what is 
the optimum experiment? What is the ideal angular resolution and survey size? 
What instrumental sensitivity is required? This article will address these questions 
(although fall a bit short of providing a complete answer). 

If we are interested only in the gravitational-wave signature, we can focus on 
the model-independent curl component of the polarization produced ho^ gravita- 
tional waves. In this article, we summarize work reported in Jaffe et al.,E3 a paper 
that extends the work of Kamionkowski and KosowskyEa and Lesgourgues et al.c2l.§ 
We ask, what is the smallest amplitude of a curl component from an inflationary 
gravitational-wave background that could be distinguished from the null hypothe- 
sis of no curl component by an experiment that maps the polarization over some 
fraction of the sky with a given angular resolution and instrumental noise? If an 
experiment concentrates on a smaller region of sky, then several things happen that 
affect the sensitivity: (1) information from modes with I < 180/6* (where 9 2 is the 
area on the sky mapped) is lost;^(2) the sample variance is increased; (3) the noise 
per pixel is decreased since more time can be spent integrating on this smaller patch 
of the sky. 

More concretely, suppose we hypothesize that there is a C component of the 
polarization with a power spectrum that has the I dependence expected from infla- 
tion, as shown in Fig. 0, but an unknown amplitude T. We can predict the size of 
the error that we will obtain from the ensemble average of the curvature of the like- 
lihood function (also known as the Fisher matrix)EKl. In such a likelihood analysis, 
the expected error on the gravitational- wave amplitude T will be <tt, where 



1 f dC ^\ 2 1 



r 2 



(9) 



§ There is also related work in KinneyEl, Zaldarriaga et al.El, and Copeland et al.0 in which it is 
determined how accurately various cosmological apdipinationary parameters can be determined in 
case of a positive detection. Magueijo and HobsorJ2£rE3 presented arguments regarding partial-sky 
coverage for temperature maps analogous to those for polarization maps presented here. 
^This is not strictly true. In principle, as usual in Fourier analysis, less sky coverage merely limits 
the independent modes one can measure to have a spacing of 81 > 180/8. In practice, instrumental 
effects (detector drifts; "1/f" noise) will render the smallest of these bins unobservable. 



8 Detection of Gravitational Waves from Inflation 



Here, the af c are the expected errors at individual J? for each Cf c multipole mo- 
ments. These are given by (cf., Kamionkowski et al.E3) 

where w^ 1 = 47rs 2 /(£ p i x iVpi x T ( 2 ) is the variance (inverse weight) per unit area on 
the sky, / s k y is the fraction of the sky observed, and i p i x is the time spent ob- 
serving each of the -/V p j X pixels. The detector sensitivity is s and the average 
sky temperature is To = 2.73 /iK (and hence the Cf c are measured in units 
that have been scaled by To). The inverse weight for a full-sky observation is 
ur 1 = 2.14 x 10 _15 t~ r 1 (s/200 /jK v / sec) 2 with t yr the total observing time in years. 
Finally Bg is the experimental beam, which for a Gaussian is Bg — e~ £ a ol 2 . We 
assume all detectors are polarized. 

The error to Cf c has two terms, one proportional to Cf c (the sample variance), 
and another proportional to w~ 1 (the noise variance). There are several compli- 
cations to note when considering these formulae: 1) We never have access to the 
actual C^ c , but only to some estimate of the spectra; 2) the expressions only deal 
approximately with the effect of partial sky coverage; and 3) the actual likelihood 
function can be considerably non-Gaussian, so the expressions above do not really 
refer to "lcr confidence limits." 

Here, we are interested in the detectability of the curl component; that is, what 
is the smallest gravitational-wave amplitude that we could confidently differentiate 
from zero? Toy problems and experience give us an approximate rule of thumb: 
the signal is detectable when it can be differentiated from the "null hypothesis" of 
Cf = 0. Thus, the £ component of the gravitational-wave signal is detectable if 
its amplitude is greater than 

af c = V2/(2^+l)/ s V>"V 2 < (11) 

We then estimate the smallest gravitational-wave amplitude T that can be distin- 
guished from zero (at "1 sigma") by using Eq. (||) with the null hypothesis Cf c = 0. 
Putting it all together, the smallest detectable gravitational-wave amplitude (scaled 
by the largest consistent with COBE) is 

where 

(2^+l)(Q cc )V 2f2 < (13) 

i>(180/8) 

Results of the calculation are shown in Fig. || Plotted there is the smallest 
gravitational-wave amplitude T detectable at 3cr by an experiment with a de- 
tector sensitivity s — 10 /iK-^sec that maps a square region of the sky over a 
year with a given beamwidth. The horizontal line shows the upper limit to the 
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width of survey (degrees) 

Fig. 3. The smallest gravitational-wave (tensor) amplitude T that could be detected at 3cr with 
an experiment with a detector sensitivity s = 10 (iKy/scc that runs for one year and maps a square 
region of the sky of a given width. The result scales with the square of the detector sensitivity 
and inversely with the duration of the experiment. The curves are (from top to bottom) for fwhm 
beamwidths of 1, 0.5, 0.3, 0.2, and 0.1 degrees. The horizontal line shows the upper limit to the 
gravitational-wave amplitude from COBE. 
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Fig. 4. Regions in the r-n parameter space space occupied by various inflationary models, as 
well as those regions that could be detected by various CMB experiments. Here r measures the 
gravitational-wave amplitude (or alternatively, the energy-scale of inflation) and n the spectral 
index for primordial perturbations. Adapted from Kinneji22l. See text for more details. 



gravitational-wave amplitude from COBE. The curves are (from top to bottom) for 
fwhm beamwidths of 1, 0.5, 0.3, 0.2, 0.1, and 0.05 degrees. The results scale with the 
square of the detector sensitivity and inversely to the duration of the experiment. 

The sensitivity to the gravitational-wave signal is a little better with an 0.5- 
degree beam than with a Tdegree beam, but even better angular resolution does 
not improve the sensitivity much. And with a resolution of 0.5 degrees or better, 
the best survey size for detecting this gravitational-wave signal is about 3 to 5 
degrees. If such a fraction of the sky is surveyed, the sensitivity to a gravitational- 
wave signal (rms) will be about 30 times better than with a full-sky survey with 
the same detector sensitivity and duration (and thus 30 times better than indicated 
in Kamionkowski and KosowskyE3E3. Thus, a balloon experiment with the same 
detector sensitivity as MAP could in principle detect the same gravitational-wave 
amplitude in a few weeks that MAP would in a year. (A width of 200 degrees 
corresponds to full-sky coverage.) 



Detection of Gravitational Waves from Inflation 11 



Since the gravitational-wave amplitude is related to the energy scale of inflation, 
Fig. |^ determines the inflationary energy scale accessible with any given experiment. 
Some indication of the range of inflationary models that can be probed with past, 
current, and future experiments is provided in Fig. ^, adapted from KinneyO. 
The parameter r (y axis) increases with increasing gravitational-wave amplitude, 
or alternatively, with the energy scale of inflation. The shaded regions show the 
where the predictions for various classes of inflationary models (e.g., exponential, 
power-law, etc; for more details see Kinney et al. The scored region labeled 

"COSMIC VARIANCE (no pol)" shows the region of parameter space that would 
be consistent with a null search for gravitational waves without polarization, while 
that labeled "PLANCK (with pol)" shows regions of parameter space that would 
be consistent with a null search for the curl component in the Planck satelliteO, an 
ESA CMB mission scheduled for launch in 2007 (in both cases, it is assumed that 
n = 0.95). 

How much could the sensitivity be improved with a post-Planck dedicated po- 
larization experiment? Achieving the Planck sensitivity will be no small feat for 
experiment, and improvements will require considerable ingenuity. Still, there are 
prospects for improvements. Very conservatively, a factor-of-3 improvement over 
Planck's detector sensitivity is plausible. The dark ellipse labeled "0f w hm = 10°, 
Cpixci = is the error ellipse that could be obtained by a putative experiment 

with roughly a factor-of-3 improvement to Planck's detector sensitivity, assuming 
that the true gravitational-wave amplitude and spectral index lie at the center of 
that ellipse. 

There are good reasons to believe that technological developments in the next 
few years may allow further improvements in detector sensitivity, perhaps of an 
order of magnitude over that achieved in Planck. As an example, we mention 
LAMB (Large-Format Array of Microwave Bolometers )E3, a new detector concept 
that would allow roughly an order-of-magnitude improvement over Planck's detector 
sensitivity with a much smaller instrument. If we assume that such a detector can be 
developed and flown in an all-sky survey, the factor-of-ten improvement would allow 
us to access the regions of parameter space that lie above the line labeled "LAMB" 
in Fig. ^. If this experiment additionally spent its time surveying a smaller region 
of the sky, then the regions of inflationary parameter space that could be accessed 
would be those that lie above the line labeled "LAMB plus scan strategy." 

5. Conclusions 

We have carried out calculations that will help assess the prospects for detection 
of the curl component of the polarization with various experiments. Our results 
can be used to forecast the signal-to-noise for the gravitational-wave signal in an 
experiment of given sky coverage, angular resolution, and instrumental noise. Of 
course, the "theoretical" factors considered here must be weighed in tandem with 
those that involve foreground subtraction and experimental logistics in the design 
or evaluation of any particular experiment. These usually encourage increasing the 
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signal-to-noise and sky coverage to better isolate experimental systematics. 

In contrast to temperature anisotropies which show power on all scales [i.e., 
£(£ + l)Cg ~ const], the polarization power peaks strongly at higher i. Hence the 
signal-to-noise in a polarization experiment of fixed flight time and instrumental 
sensitivity may be improved by surveying a smaller region of sky, unlike the case 
for temperature-anisotropy experiments. The ideal survey for detecting the curl 
component from gravitational waves is of order 2-5 degrees, and the sensitivity is 
not improved much for angular resolutions smaller than 0.2 degrees. An experiment 
with this ideal sky coverage could improve on the sensitivity to gravitational waves 
of a full-sky experiment by roughly a factor of 30. When coupled with realistic 
forecasts for improvements in detector sensitivity, we find that an experiment that 
accesses a very good fraction of the inflationary parameter space (specifically, most 
of the inflationary parameter space associated with grand unification) is conceivable 
in the not-too-distant future. 

Before closing, we should note that secondarv_(in the density-perturbation am- 
plitude) effects suctLas weak gravitational lensingcj or re-scattering of CMB photons 
from reionized gasCJ may lead to the production of a curl component in the CMB, 
even in the absence of gravitational waves. However, these secondary effects should 
be distinguishable from those of gravitational waves, as they produce a curl compo- 
nent primarily at angular scales much smaller than those at which the gravitational- 
wave signal should show up. Of course, an angular resolution better than that we 
have suggested here and a survey area a bit larger than we have suggested here may 
be required to distinguish the gravitational-wave signal from these other sources of 
a curl component (as well as from foregrounds) . A more complete assessment of the 
impact of these secondary effects on the detectability of gravitational waves is now 
underwayo. 

Finally, the CMB polarization will be useful for a wide variety of other purposes 
in cosmology. For example, detection, and ultimately mapping, of the polariza- 
tion will help isolate the peculiar velocity jit the surface of last scattered, con- 
strain the io nizati on history of the Universec3, deter rnine th e nature of primordial 
perturbationsOEJ, probe primordial magnetic fieldsOEHcj and cosmologicaJ par- 
ity violationEjEH, and maybe more (see, e.g., Kamionkowski and KosowskytJ for a 
recent review). 
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